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ABSTRACT 

The hard X-ray (> 2 keV) emission of the local and distant Universe as 
observed with the HEA01-A2 experiment is reconsidered in the context of large 
scale cosmic structure. 

Using all-sky X-ray samples of AGN and galaxy clusters we remove the 
dominant local X-ray flux from within a redshift of ~ 0.02. We evaluate 
the dipolar and higher order harmonic structure in 4 X-ray colours. The 
estimated dipole anisotropy of the unresolved flux appears to be consistent 
with a combination of the Compton- Getting effect due to the Local Group 
motion (dipole amphtude A = 0.0042) and remaining large scale structure 
(0.0023 ^ A <, 0.0085), in good agreement with the expectations of Cold Dark 
Matter models. 

The observed anisotropy does however also suggest a non-negligible Galactic 
contribution which is more complex than current, simple models of > 2keV 
Galactic X-ray emission. Comparison of the soft and hard colour maps with a 
harmonic analysis of the 1.5 keV ROSAT all-sky data qualitatively suggests that 
at least a third of the faint, unresolved ~ 18° scale structure in the HEA01-A2 
data may be Galactic in origin. However, the effect on measured flux dipoles is 
small ( < 3%). 

We derive an expression for dipole anisotropy and acceleration and 
demonstrate how the dipole anisotropy of the distant X-ray frame can constrain 
the amplitude of bulk motions of the universe. From observed bulk motions over 
a local ~ 50 h^^ Mpc radius volume we determine 0.14 ^ QQ'^/bx{0) ^ 0.59, 
where Qq is the universal density parameter and 6a; (0) is the present epoch bias 
parameter, defined as the ratio of fluctuations in the X-ray source density and 
the mass density. 

Subject headings: X-rays: general — cosmology: observations, large-scale 
structure of universe 



- 3 - 



1. Introduction 



Establishing all the sources of the extragalactic flux dominating the X-ray sky remains 
a fundamental challenge. In the soft X-ray band (0.5 — 2 keV) some 70-80% of the observed 
flux can be accounted for by extrapolation of those objects resolved in deep fields ( pasinger 



et al 19981) . Such objects are QSO's, AGN and possibly narrow emission line galaxies. In 
the hard X-ray band ( ^ 2 keV), however, the situation is less clear. Deep surveys with 
ASCA, reaching 1 x 10~^^ergs~^ cm~^ in the 2-10 keV band now account for ~ 30% of the 
XRB ( [Cagnoni et al 19^ , |Ueda et al 1999|) . Although recent results with Chandra (e.g. 
Mushotzky et al 2000| ) have now pushed this to ~ 75%. In addition, the spectral form of 
the still unresolved flux does not fit with the spectra of any single class of known objects, 
although more detailed models have had some success ( [Leiter fc Boldt 1992| , [Comastri et| 
al 1995| , pVladau et al 199^ ). The investigation of X-ray emission associated with the local 
Universe has produced more tangible results. The autocorrelation function ( [Jahoda 1993| ) 
of the unresolved X-ray flux and its cross-correlation with other extragalactic catalogues 
( [Jahoda et al 199"T| , [Jahoda et al 199^ , [Lahav et al 1993| , [Miyaji et al 1994[ , Parcons et al 
1995[ [Refregier et al 1997[ , [Almaini et al 19971) provides useful information on the volume 



emissivity of that fraction of X-ray emission correlated with galaxies. 

For such studies the HEA01-A2 all-sky survey continues to be the best hard band 
all-sky data available. The integrated flux in this survey originates from 2; = to a 
redshift of 2; ~ 4; and might therefore be considered the most complete survey of baryonic 
matter currently available at any wavelength. Less than 2% of the total extragalactic flux 
comes from sources identified in existing all sky catalogs ( piccinotti et al 1982 , Grossan 
1992[) . Since the majority of the observed flux comes from high- 2; we expect it to be 



highly isotropic, at least to C(10~^), with possible deviations caused by anisotropies in the 
population of nearby but unresolved sources. Indeed, the extragalactic hard X-ray emission 
associated with foreground sources within z ~ 0.02 is highly anisotropic, indicative of the 
pronounced structure in the mass distribution traced by present-epoch AGNs. For example, 
the dipole to monopole ratio observed in this population ( [Miyaji fc Boldt 1990[ , [Miyaji 
et al 19911 , [Miyaji et alTUd^) is very large, ~ 0.5, while the monopole contribution to the 
total extragalactic emission is ~ 1%. The observed AGN dipole has also been demonstrated 
in these studies to be compatible with the direction and magnitude of the local group 
velocity. This is a crucial observational 'calibration' which supports the idea that AGN 
X-ray sources do indeed trace the underlying gravitational mass distribution responsible for 
peculiar motion. 

Lahav, Piran & Treyer (1997) (hereafter LPT97) have evaluated the expected 
large-scale angular fluctuations in the XRB for a range of power spectra of mass fluctuations 
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(e.g. Cold Dark Matter (CDM) models) and X-ray evolution scenarios. In a foUowup 
paper, and companion to this present work ( [Ireyer et al 199^ ) we have investigated these 
large-scale (10° — 180°) fluctuations. Using the HEA01-A2 data the power spectrum of mass 
fluctuations can be probed on scales of ~ 600h~^ Mpc (where h is the present epoch Hubble 
constant in units of 100 km s~^ Mpc~^). This data can also be used to constrain the fractal 
correlation dimension {D2) of structure on this scale ( [Peebles 1993| , [Wu et al 1999|) and 



hence test the homogeneity of the Universe and the validity of the Cosmological Principle. 
A value of -D2 = 3 to a precision of 10~^ is observed, strongly supporting homogeneity on 
large (~ 600h~^Mpc) scales. 



Work has also been done ( [Shafer 1983| , [Boldt 19E1\ , [Jahoda et al 1992] , [Jahoda 1993D to 



determine the extragalactic X-ray flux dipole from the HEAO data. The motivation for this 
being the determination of the expected Compton-Getting (CG) dipole due to our motion 
with respect to the distant X-ray frame, which is expected to agree with the direction and 
velocity inferred from the CMB dipole ( [Lineweaver et al 1996|) . The amplitude of such 



a dipole constrains the distance of the frame which can be considered to be at rest with 
respect to the CMB. 

However, as demonstrated in LPT97, for a typical observer in (for example) a CDM 
universe the amplitude of the expected CG dipole and that due to emission correlated with 
large-scale structure (LSS) are comparable. This coupling of the two dipole terms makes 
it hard to use the XRB dipole to conflrm the CMB motion dipole. Other investigations 
( [Plionis fc Georgantopoulos 1999|) of the XRB dipole using the soft, l.SkeV ROSAT 
all-sky survey data seem to provide some conflrmation of this prediction. However, as we 
demonstrate below, data in this softer band is strongly contaminated by Galactic emission 
compared to the hard, > 2 keV HEAO data. 

In this present work we reconsider the HEA01-A2 data and examine the X-ray 
background dipole together with higher order structures and investigate the relative 
importance of the CG effect, LSS, Galactic emission and the relationship of dipolar 
structure to bulk motions. In Section 2 we describe the HEA01-A2 all sky data, and a 
newly discovered small instrumental effect and a prescription for its removal. In Section 3 
we describe and apply spherical harmonic analysis to the large angular scale structures in 
the data and investigate the effects of removing known sources. In Section 4 we present 
angular power spectra. In Section 5 we measure the flux dipole of the unresolved XRB 
and assess its signiflcance through simple simulations. In Section 6 we derive the full 
cosmological expression for dipole anisotropy and its relationship to peculiar velocity. In 
Section 7 we discuss observations of local bulk motions and apply these results to our X-ray 
dipole estimates to obtain constraints on bias parameters. In Section 8 we summarise our 
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results and present conclusions. 

2. The HEA01-A2 data 

We have taken the present data from the online archives at the High Energy 
Astrophysics Science Archive Research Center (HEASARC) at the NASA/Goddard Space 
Flight Center Q. In its raw state the data used here consists of the 'Small Field of View' 
(SFOV) surface brightness in counts/sec per beam, stored in 0.25° x 0.5° (a total of 
720 X 720) pixels in a rectangular projection in ecliptic coordinates. The intrinsic resolution 
size of independent data points is however 1.5° x 3°. In addition, for the purposes of 
estimating the instrumental background, we have utilized the equivalent 'Large Field of 
View' (LFOV) data, which although stored in similar format (720 x 720 pixels) has an 
intrinsic resolution of 3° x 3°. The all-sky survey is available in 4 overlapping energy 
bands or colours: Soft, Hard, Total and R15. Allen et al (1994) presents effective area 
curves (detection efficiency as a function of energy) for these bands. The Soft color consists 
primarily of photons detected in the first layer of the argon and xenon filled detectors and 
with pulse height less than 6 keV; the Hard color consists of second layer photons and large 
pulse heights from the first layer and has very little response below 6 keV while the soft 
color has very little response above 8 keV. The R15 color was, throughout the mission, the 
unweighted sum of x-rays detected in the first and second layer of the High- Energy-Detector 
(HED)-3 and the second layer of the Medium-Energy-Dectector. The weights which define 
these colours are chosen so that a source characterized by a photon index of -1.7 which 
produces 1 R15 count/sec per beam will also produce 1 Total (or 1 Soft or 1 Hard) count/sec 
per beam. The R15 band is the most stable colour over the period of observations, the 
Total band has the highest signal-to-noise. The 4 bands have effective areas peaking at 3, 
7, 6 and 10 keV for Soft, Hard, Total and R15 respectively. All data was taken from the 6 
month observation period beginning day 322 of 1977. To convert the raw surface brightness 
data into standard units we use a conversion factor of 2.2 x 10~^^ erg s~^ cm~^ (4.5 deg^)~^ 
per count/sec. The beam area (4.5 deg^) reflects the instrinsic 1.5° x 3.0° resolution of the 
SFOV. 

In Figure 1 the raw HEA01-A2 data is presented in a Hammer-Aitoff, Galactic 
projection. The data consists of one complete scan of the sky using the Total band ( [Alien 
let al 19941 ). 

The darkest pixels are associated with individual bright sources and diffuse Galactic 



{ ftp://legacy.gsfc.nasa.gov/heaol/data/a2/maps/heasarc_med_hed/ 



) 
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emission. Only ~ 100 of the brightest sources at high galactic latitudes can be identified 
( [Piccinotti et al 1982| ) identifies 17 galactic and 68 extragalactic sources). In all of the 
following we restrict ourselves to |6| > 20° where complete catalogs are available. 

The data is first corrected for the locally estimated instrumental background. This is 
achieved by using both the LFOV and SFOV data to solve the following two equations at 

each pixel for Ibackground- -^sfov -^background ~l~ -^sky and Iljov -^background ~l~ wherC 

the factor of 2.26 is the ratio of area solid-angle products for the LFOV and SFOV. To 
reduce noise in the background estimate we then correct the data by subtracting the mean 
background estimated in strips of constant ecliptic longitude. Prior to the background 
calculation we remove sources and mask Galactic regions as described below. 

In the course of examining the data we observed a systematic change in the measured 
fiux as a function of observation date. Since the data was taken in great circles through 
the ecliptic pole, data separated by 180° in ecliptic longitude correspond to approximately 
the same epoch of observation. Taking the mean fiux at fixed longitude of the source 
removed/masked data (see below), wrapped by 180° we discovered a clear linear trend 
in the observed counts. In Figure 2 we plot the observed effect in the Total band 
data. The least squares linear models for all 4 HEAO colours are shown in Figure 
3. The slopes of the trends (in units of counts/sec, vs longitudinal pixel index) are 
-0.000716, -0.000194, -0.000357, -0.000265 for the Soft, Hard, Total and R15 colours 
respectively, with maximum and minimum at longitudinal pixel indices 69 and 429 
respectively. 

The time dependent term therefore ranges from — 7% of the mean intensity for the 
Soft band, and less for the other colours. In all subsequent analysis the HEAO data has 
been corrected by subtracting the least squares trend. The resulting change on all quantities 
discussed is small and always less than 10%. 

To improve computational efficiency in all the following analysis we have re-binned the 
data (following background, and systematics correction, and source removal as described 
below) into 3° x 6° pixels in ecliptic coordinates (smaller resolution pixels are more strongly 
correlated due to the instrument beam) 

3. Large scale structures in the HEAO data 

Unlike many all-sky catalogues (e.g. IRAS, optical surveys etc.) the ability 
to unambiguously separate foreground (Galactic/local emission) from background 
(extragalactic) information is limited in the HEAO X-ray data. The total number of 
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resolved foreground and background sources is small (~ 100) and a detailed model of 
possible large scale Galactic emission is hard to determine. The Galactic 2 — 60 keV 
emission model of Iwan et al (1982) predicts that some 2% of the observed emission at the 
galactic poles in the A2 band is of Galactic origin. The largest contribution at latitudes 
b > 20° is 5%. (The Iwan model predicts galactic contributions of up to 10% at low 
latitudes, although there is certainly another, more centrally concentrated component as 
well ( [Worrall et al 1982| , [Warwick et al 1985| , [Valinia fc Marshall 199"^ ). More recently, 
studies in the soft bands (< 0.75 keV) by ROSAT ( pnowden et al 19971) indicate that, at 
these lower energies, the picture is more complicated, with structure at all scales. Whether 
this soft emission distribution is a good indicator of the much harder > 2 keV emission 
is unclear; the galactic contribution in this soft band is almost certainly larger than that 
above 2 keV, and is potentially more variable as well. In this present work we attempt 
to at least qualitatively evaluate the likely foreground vs. background contributions to 
further constrain our estimates of the extragalactic flux anisotropy. In order to evaluate the 
structure in the map of Figure 1 we use spherical harmonic analysis to filter the high (noisy) 
frequencies from the map and to reconstruct the flux variations (e.g. Scharf et al 1992). 
Briefly, the map is expanded into the orthonormal set of spherical harmonics {Yim,{0 , (p)) by 
determining the harmonic coefficients as a sum over the flux cells: 



aim = Y.IiAuiY^{ui) , (1) 

i 

where /j is the mean surface brightness in a cell, Aui is the cell area in direction Ui. The 
surface brightness at any point can then be reconstructed using only lower order harmonics, 
where the resolution is determined by the highest harmonic and goes as ~ Tr/lmax- In our 
definitions below the monopole Mx is defined as Anl thus aoo = V^ttI = M^/V^- 



3.1. Higher order anisotropies 

In Figure 4 a harmonic reconstruction of the raw (Total band) HEAO data is shown 
to a resolution of ~ 18° {Imax = 10). The data is clearly dominated by emission associated 
with the Galaxy (either resolved or unresolved sources). As a first step towards removing 
this foreground we construct a 'mask' using the list of resolved and identified Galactic 
X-ray sources ( piccinotti et al 1982| ) and a |6| < 22° Galactic Plane mask. Regions of sizes 



varying from ~ 8° diameter to 12° diameter are excised around resolved sources, larger 
regions are removed around the Large and Small Magellanic Clouds (LMC, SMC). A total 
of ~ 38% of the sky area is removed by this mask. The results are dramatic, in Figure 5 
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the harmonic reconstruction (to l^ax — 10) of the data is shown in all four bands, with the 

above Galactic masking applied. Those regions excised have been filled with the new mean 
flux as a first order correction (see discussion of dipole estimation below). As described in 
Section 2, the count/sec in each colour are weighted to be equivalent for emission with a 
photon index of -1.7. The observed differences in the structure between (for example) the 
Soft and Hard bands can then be directly interpreted as spatial differences in the mean 
spectral index of X-ray emission (modulo variations in the signal-to-noise) . 

Two strong flux enhancements are apparent. The uppermost (at b ^ 80° and spanning 
the Northern Galactic cap) we associate with the Virgo and Coma clusters. The lower peak 
(at / ~ 315° and b ~ 30°) is close to the Centaurus/Great Attractor region (e.g. Scharf et 
al 1992, Webster et al 1998). However, we note that it is also closer to the Galactic Plane 
and may not be free of Galactic contamination. 

Next we make use of the Pincinotti (1982) catalogue of extragalactic sources 
(supplemented with the clusters of Edge et al (1990) with z < 0.003) and excise these 
regions (in ~ 8° diameter cuts) in addition to the above Galactic masking. This results in 
a dataset with all sources removed down to a flux limit of 3 x 10~^^ ergs s~^ cm~^ (c.f. 
Treyer et al 1998). The median redshift of sources in the extragalactic Piccinotti sample 
is known to be ^ ~ 0.02. In a further effort to remove all signiflcant X-ray sources from 
within the local volume, we have made use of the HEAO Al catalogue of source detections 
(Ron Remillard, private communication). Using this data we have removed all sources to 
a slightly lower flux limit of 2 x 10^^^ ergs s~^ cm~^ (2-10 keV). We determine this flux 
limit value from both our chosen limiting count rate in the Al catalogue (0.006 ct/s) and 
by normalising with respect to the extrapolated Piccinotti LogN-LogS. This yields a final 
monopole value corresponding to ~ 6.2 x 10^^ erg s ~^ Mpc~^. The final unmasked 
sky area is 48%, in Figure 6 the corresponding harmonic reconstructions (to l^ax = 10) are 
shown. 

Prom the X-ray luminosity functions of Grossan (1992) and Boyle et al. (1998), we 
estimate that the mean luminosity of the local X-ray source population is ~ 5 x 10^^ erg 
s~^, assuming a lower luminosity cutoff of 10^^ erg s"^ and a present day Hubble constant 
Ho=100 kms~^Mpc~^. As the luminosity function is dominated by faint objects and and 
the low-luminosity cutoff is somewhat arbitrary, is likely to be a lower limit. Removing 
sources brighter than 2 x 10^^^ erg s~^ cm~^ corresponds to removing all sources brighter 
than out to 2; ~ 0.015, and all sources brighter than 10 x out to 2; ~ 0.05. Therefore 
most X-ray sources with z > 0.015 are still unresolved and contribute to the measured 
dipole (see below). 

The contour levels in Figure 6 are approximately a factor ~ 3 more exaggerated relative 
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to the mean than in Figure 5. Unhke the extragalactic sky dominated by resolved sources, 
significant differences in structure is now seen between the Soft and Hard colours. Notably, 
in the Soft band two strong patches are seen at (~ 10°, ~ —28°) and (~ 80°, ~ —40°), 
and a significant structure at / ~ 180°, extending over 15° ^ b ^ 80°. In the Hard band 
the most significant features are seen at (~ 200°, ~ —40°) and (~ 280°, ~ 60°). Possible 
identifications with known structures are given in the Figure 6 caption. 

In Figure 7 we plot the harmonic reconstructions of 'noise' skies for both the Galactic 
and full extragalactic mask cases in Figures 5 & 6, to provide a qualitative guide to the 
significance of features. From these figures wc can assess the signal-to-noise of structures 
seen in Figures 5 & 6 (Total band). The most significant individual features are typically 
a factor 2-3 higher in counts/sec when compared to the equivalent peaks in the pure noise 
levels. It is apparent that there are at least two likely real, positive, features in Figure 5 
and possibly four in Figure 6. 

In an attempt to assess the likely amplitude and morphology of Galactic contamination 
in the HEA01-A2 data we have analysed the all-sky X-ray data from the ROSAT All Sky 
Survey (RASS) using the maps of Snowden et al (1997, 1995). In Figure 8 we present 
a harmonic reconstruction of the 1.5 keV RASS data, subject to the same Galactic and 
extragalactic source mask as the HEA01-A2 data as shown in Figure 6. The RASS data 
was used in a 480 x 240 pixel Hammer-Aitoff Galactic projection as available from the 
online archives, no further processing was made, other than filling empty pixels with the 
mean fiux after applying the HEA01-A2 Galactic and extragalactic source masks. As 
discussed in Snowden et al (1995) the 1.5 keV band data spans the energy range ~ 0.7 — 2 
keV, with essentially no counts above 2 keV. The HEA01-A2 Hard band is therefore the 
only colour with no response in this energy range. The response in the other bands is also 
modest at 2 keV, the effective area is down from its maximum by factors of ~ 3, 10, and 
10 for the Soft, Total, and R15 colours. The three most prominent features in the RASS 
reconstruction are concidcnt with the features identified in Snowden et al (1995) as the 
Northern Spur/Loop 1 (/ ~ 20°, b ~ 40°), and the northern and southern parts of a Galactic 
Bulge emission (/ - 340°, b - 30°, / - 0°, 6 30°). 

The location of the two Bulge structures is particularly worrying for our interpretation 
of the structures in Figure 6. In particular the southern structure seen most strongly in 
the Soft band at / ~ 10°, b ~ —30° is not present in the Hard band data and might indeed 
be correlated with the RASS southern Bulge feature. Snowden et al (1997) have further 
investigated the RASS maps and propose that a soft Galactic bulge component is well 
modelled by a cylindrical volume (with exponential density fall-off away from the Galactic 
plane) of ~ gas. As discussed above, this does not preclude Galactic bulge emission 
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above 2 keV, but for a thermal spectrum very little emission would be expected. The 
southern HEAO 'Bulge' feature appears to fall into this category. The HEAO structure 
at / ~ 325°, 6 ~ 30° is however present in all 4 colours, which suggests that while it could 
be contaminated by the northern Bulge feature, it is more likely to have an extragalactic 
origin. Therefore, we conservatively estimate that at least 1 out of 3 of the dominant faint 
(soft), unresolved structures in the HEA01-A2 data may be Galactic in origin. This is 
consistent with the detailed investigation of the RASS 1.5 keV dipole made by Plionis & 
Georgantopoulos (1999) who estimated that 20 — 30% of the total background is Galactic 
in origin in that data. 

Without much more rigorous cross-comparison any correlations of the RASS with 
HEA01-A2 should be made with caution, since the intrinsic spatial resolution of the RASS 
(~ 2°) is significantly different to that of the HEA01-A2. In Section 5 below we investigate 
the effect on the flux dipole of Galactic emission as described by the Iwan et al (1982) 
model, as well as the effect of removing the putative 'Southern Bulge' feature. 



4. Angular power spectra 

In an associated paper ( P?reyer et al 1998|) we have investigated the angular fluctuations 
of the XRB in the context of models of large scale structure (as formulated by LPT97). 

Following the formalism used in those papers the total predicted angular power (taken 
as an ensemble average) is a combination of a large-scale structure component (signal) and 
a component due to the discreteness of sources (shot-noise) : 



(|a-|n =(|arr) +(\aT?) (2) 

\' ' ' /model \l ' I / LSS V '■ ' / SN ^ ^ 

This is then convolved with the foreground mask to produce the observed angular 
power spectrum: 



{\cT\') = T.\wr'\'{wf\') (3) 

I'm' 

Finally we normalize over the monopole, and for total signal use the notation: 



i{\cr\')Lss+{\cr?)sNy^' 



(l«ooP) 



1/2 



(4) 



- 11 - 



(see Section 3). 

Figure 9 plots the normalized angular power spectra (i.e. relative to the monopole, 
Treyer et al 199^ , Eqn 14) of the HEAO data for the complete extragalactic plus Galactic 
mask in all 4 colours. The expected CG effect (see §5 below) has also been subtracted from 
the data. (Note that in Treyer et al (1998) fewer sources were removed, corresponding to 
a slightly higher flux limit of 3 x 10^^^ erg s~^ cm~^ (2-10 keV)). The Total band spectra 
used in Treyer et al (1998) are consistent with those shown here, despite some differences in 
data processing and source removal. The solid curve in Figure 9 is the predicted Total band 
angular spectra for a low density CDM model with pure X-ray source density evolution and 
constant, linear biasing (see §6 below and Treyer et al (1998)), and including the predicted 
noise due to unresolved sources (LPT97). 

It is reassuring that the 4 spectra are in generally close agreement (despite different 
noise characteristics). We do however note that in Figure 9, the Soft band spectra (possibly 
subject to the most Galactic contamination) shows a strongly discrepant quadrupole (/ = 4) 
term. 

In Figure 10 we also plot the mean spectra obtained from 'noise' simulations. As 
described in Section 5 below these consist of sky fluxes randomly drawn from the real data 
and masked as per the real data. The complete sky mask is applied both to the data from 
which fluxes are drawn and to the simulations. The noise estimated from simulations is in 
agreement with the analytic estimate given in LPT97 and Treyer et al. (1998). 

It is apparent that the harmonic shot noise level varies between the HEAO colours. 
As expected the Total band has the lowest noise level, followed by the R15 band and the 
highest noise in the Soft and Hard bands, at very similar levels. Comparing these spectra 
with those in Figure 9 we can see that by / ~ 5 in the real data we are approaching the shot 
noise level and by / ~ 10 we have practically reached the shot noise level. 

5. The observed flux dipole 

The dipole flux anisotropy is described by the harmonic coefficients aim, and can also 
be written or derived as a simple vector D^,, where; 

= ^ liAuiTi (5) 

i 

We can of course parameterize this differently. For example, the CG effect due to our 
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motion with respect to a radiation field (c.f. CMB dipole etc) is a purely dipolar effect 
of the form /' = /(I + Acos0) where I is surface brightness, and 9 is measured from the 
direction of motion, in which case ID^-I/M^ = A/3 (note that = AttI). The factor A for 
the 2-10 keV X-ray band CG effect is: 

A = -[l + («/3)], (6) 



where the energy spectral index is a = 0.4 ( [Boldt 19"S^ ). Assuming the Solar velocity 
with respect to the CMB ( Lineweaver et al 1996 ), we determine A = 4.2 x 10~^ for the 



X-ray CG effect. Here we choose to present the dipole amplitude as A using the above 
conversion from D^. 

If an all-sky, foreground removed. X-ray surface brightness map existed the flux 
dipole could be immediately obtained as a vector sum over the flux cells. However, as 
demonstrated above, removing the foreground involves removing information about the 
background too. The first order correction to the removal (masking) of regions when the 
harmonic decomposition is performed is to fill those regions uniformly with the mean 
density/fiux of the unmasked regions. This will not, however, remove cross-talk between 
the true (full sky, orthonormal) harmonic coefficients (including the dipole) (e.g. Scharf et 
al (1992)). One solution is to attempt to reconstruct the full sky harmonics by inversion of 
the coefficient matrix, suitably controlled by (for example) a Wiener filter (e.g. Lahav et 
al (1994)). While this is a very powerful method it does require a model of the expected 
harmonic power and a full knowledge of the noise matrix. Additionally, such reconstruction 
is severely limited by the amount of masking (for a realistic harmonic spectrum of large 
scale structure), in the case of the HEAO data this is large; if all resolved Galactic and 
extragalactic sources and the Galactic plane are masked ~ 52% of the sky is removed. 

In this present work we use two different methods of dipole estimation. First, we 
perform the vector sum over all flux cells, with a first order correction of filling masked cells 
with the mean flux over the unmasked regions (the spherical harmonic approach, hereafter 
Method 1). Second, we perform a least squares fit of a cos{9) dependence dipole to the 
data and determine the best fit {l,b) and A (hereafter Method 2). This latter method 
does not make use of the masked region, and does not suffer from cross-talk, but it does 
assume the specific form of the dipole anisotropy, and (unrealistically) that the residuals 
are negligible. In the following we apply it only to the Total band data, since this data has 
the best signal-to-noise, and evaluate both methods using Monte-Carlo simulations below. 



The results of the various dipole estimations for the HEAO data are presented in 
Tables 1 and 2 for Methods 1 and 2 respectively. 
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As a comparison with previous works on the dipole moment of X-ray AGN ( |Miyaji fc 



Boldt 1990| , [Miyaji et al 1991| , [Miyaji 1994| ) we subtract the Total band, source removed 



extragalactic dipole (row 7, Table 1) from the full extragalactic Total band dipole (row 3, 
Table 1). The resulting vector (which is equivalent to the dipole vector of the removed 
extragalactic sources) has A ~ 8.2 x 10~^ and points at (316°, 53°). Miyaji (1994) found 
that the flux dipole of resolved AGN (which has different shot noise) in the HEAOl survey 
to z ^ 0.015 has a direction (293°, 33°) ± 20°, consistent with this difference vector. In their 
analysis of the RASS 1.5 keV XRB dipole Plionis & Georgantopoulos (1999) determine a 
best dipole direction of (288°, 25°) and amplitude A = 0.051 (converting their D^/M^). 
While the direction is in general agreement with our results their amplitude is more than 
a factor of ~ 10 higher compared to our results in Tables 1 & 2. This discrepancy is likely 
due to the increased difficulty of diffuse foreground removal in the soft ROSAT band (< 2 
keV) and an increased soft-band contribution from galaxy clusters and groups. Indeed, 
Plionis & Georgantopoulos (1999) estimate that Virgo contributes as much as 20% to the 
RASS dipole amplitude. 

As a demonstration of the effect of applying a correction assuming the Iwan et al 
(1982) Galactic emission model the 3rd block of Table 1 presents the results of removing a 
Iwan Galactic emission model (normalized to 3% of the monopole) from the distant flux 
dataset. The effect is relatively modest and systematically decreases the dipole amplitude 
A, and moves the observed dipole direction away from the Galactic center. We also assess 
the effect of removing the expected CG effect from the data (block 4, Table 1). The dipole 
is significantly altered; A drops by 15 — 20% and the direction changes by as much as ~ 50°, 
towards the Galactic center. Combined with a 3% Iwan et al. correction (block 5 Table 1) 
A is further reduced. We also measure the variation in A (and direction, see Figure 13) as a 
function of increasing Galactic correction (Figure 11), for cases with and without removal of 
the expected CG effect (lower and upper curves respectively). In the case in which the CG 
effect is removed (after the Galactic correction) the amplitude of the Galactic normalization 
is the most pronounced, reducing A by a factor of 4 between a 1% and 13% correction. 

In Table 2 the equivalent results from the Method 2 dipole measurements (using the 
Total band) are presented. While the directions agree fairly well with those of the Method 
1 dipole estimates, the amplitudes appear to be systematically larger by factors ~ 2.5 — 3 
(see below). 

Finally, we test the effect of removing a 30° x 30° patch of sky around the putative 
'Southern Galactic Bulge' region, centred on / = 0°,6 = —25°. As done above this region 
is filled with the mean flux. The effect on both Method 1 and 2 dipole estimates is small, 
reducing A by 3 — 5% and altering the direction by 56' ~ 10°. 
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In order to test the significance of the dipole measurements and the abihty of the 
two methods in recovering a genuine dipole signal we use simple Monte-Carlo simulations. 
Taking the fluxes of a dataset with the Galactic and extragalactic mask applied we 
resample the flux distribution and construct a random sky map (equivalent to assuming no 
correlation between flux cells) which is subject to the same sky incompleteness (in this case 
the complete mask). The results of harmonic analyses on the simulations over all bands 
and to Z = 20 have been shown previously in Figure 10. Here we concentrate on the dipolar 
measurements. 

The mean dipole amplitudes over 10 realizations, are shown in Table 3. It is clear 
that the 'noise sky' dipoles are significantly smaller than the dipoles seen in the real data, 
indicating the presence of genuine correlated structure. The bottom two rows of Table 3 
show the result of adding a reahstic CG effect to the simulated data. It is encouraging that 
the estimated amphtudes are close to the input value (A = 4.2 x 10~^), and the dipole 
directions are in general agreement, but we note that the two methods appear to differ in a 
systematic way, with the Method 1 estimate being consistently smaller by a factor 1.5 — 1.8. 
This is not unexpected. As discussed in Trcycr ct al (1998) and Scharf et al (1992) (and 
references therein) an incomplete sky creates cross-talk between the harmonic coefficients. 
In the case of the mask used here the net result is to systematically lower the observed 
amplitude of the Method 1 dipole. The Method 2 dipole estimate does not suffer from 
such an effect, although it is a less general estimate of dipole anisotropy. The difference 
seen in the simulations accounts for at least 50% of the Method 1/2 A discrepancy seen 
in the HEA01-A2 data. More detailed simulations would be needed (including reahstic 
fluctuations instead of pure noise) to determine the precise difference expected. A full 
treatment of the significance of the observed dipole is beyond the scope of this present work. 

However, on the basis of the simulations (Table 3) and the results of Tables 1 and 2, 
we estimate that our observed dipoles are significant at greater than a ~ 2 — 3(7 level, and 
have a typical direction error of ~ 30°. 

6. The flux dipole and bulk motions 

The X-ray fiux dipole observed at frequency uq is defined as: 

B, = J2fiMri, (7) 

i 

where the sum is over all directions in the sky, and /i(i'o) is the integrated X-ray fiux in the 
direction fj. 
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Following the formalism given in LPT97 and assuming linear, epoch-dependent biasing 
bx{z), such that fluctuations in X-ray sources and mass are related by 5x{rc,r) — bxS{rc,r), 
then: 

= 1 1 0(L., z)^^^^il±^[l + b,{z)S{r,, f)]f dVJL, , (8) 

where is the radial probability of a source with luminosity Li, at redshift z and S is 
the mass density contrast. 

If the number density of the X-ray sources evolves as {1 + z)''', their luminosity as 
{1 + zY, and Li, oc 1/"°', then we can define q — d + e — a + 1 and the X-ray volume 
emissivity as: 

p,{z) = I L,cf>{L,,z){l + z)dL,^p,{0){l + zy . (9) 

The dipole can then be written in the form of a "dipolar Olbers' integral" : 

J_ r pAz)b^iz)S{r,,r) 
^ At: J r^,{l + zy "^^^ ^'"^ 

Recall that in an Einstein-de Sitter Universe {VLq = 1) dVc — r^dvcduj and 
Tc — '^'r'H\\ — (1 + 2;)"^^^], where th = cj Hq is the Hubble radius. 

As wc do not have a model for 5{v) in our neighbourhood, we can only make statistical 
predictions using a model for the power-spectrum P{k) (LPT97, Treyer et al 1998). Of 
course, what wc observe is a single realization and this one realization may not be well 
represented by the rms value. The rms dipole {I = 1) can be expressed as (see Eq. 7 in 
Treyer et al 1998): 

{Wm?)LSS = j k'P{k)\^m'^K (11) 

where the window function contains the various model parameters: 

^i{k) = aMz){l + zy-^'^3i{kr,)W,ut{z)dz , (12) 

and the function Wcutiz) accounts for the removal of sources brighter than a given flux 
cutoff Scut- CTs is the usual rms mass fluctuation in an 8h~^Mpc radius sphere. As in Figure 
8, we use a fiducial model assuming low density CDM, pure density evolution with q — 4.6 
and Zmax = 1-3 (based on Hasinger 1998), and constant biasing. 

Figure 12 shows the growth of the rms flux dipole as a function of the outer radius 
cutoff Zjnax for our fiducial cosmological model and 3 scenarios for Wcut{z)- The figure 
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shows, first, that a fiux cutoff of 2 x 10"^^ erg s~^cm~^, as used in the present data analysis, 
is very similar to removing all sources within 2; < 0.01 — 0.015; and secondly, that the rms 
dipole converges very rapidly, so that most of it originates from z < 0.05 — 0.1. Consequently 
there is very little signal due to structure lying further out, and in the presence of noise we 
will have effectively no information from z > 0.1, at least in the rms sense. As the growth 
of the dipole depends on the power spectrum, P{k), in models with more large scale power 
the convergence with z will be slower. Note that this convergence is not due to an "Gibers' 
effect" . The total intensity / of the XRB keeps increasing to much higher redshift than 
do the fluctuations: / oc /cf™"(l + z)'^-'^/'^dz while ^'i(A;) oc /o''"""(l + z)i-^/'^ji{krc)dz (to 
flrst approximation). Unlike the monopole, the high redshift fluctuations (dipole and higher 
harmonics) are effectively washed out by angular averaging over the sky (governed by the 
Bessel function dependency). 

Therefore, we can only use the XRB flux dipole to constrain large scale structures out 
to 150-300 h~^ Mpc. Coupled with our estimate that the bright sources we remove are 
distributed to a distance of ~ 60 h~^Mpc, we should be able to compare our XRB dipole 
measurement with direct measurements of bulk flows over a similar volume. 

In linear perturbation theory the peculiar motion at any point in space is directly 
proportional to the gravitational acceleration, we can therefore write (assuming all motion 
was zero a Hubble time ago ( Peebles 1980|) : 

V ^ l%'-'sHo' , (13) 

where Qq is the density parameter. The gravitational acceleration g in Newtonian gravity 
is: 

g = Gp{0) J^-^idV, (14) 

(p(0) is the present-epoch mean mass density). We note that this expression only holds in 
on small scales, by choosing locally Minkowski coordinates (Peebles 1993, p. 268). We can 
therefore approximate Eqn. 10 above for the flux dipole at low redshift to be: 

= ^p.(0)6,(0)g/G' (15) 

Therefore at low redshift the flux of a source follows an inverse square law and if light 
traces mass in a spatially invariant and linear way then any anisotropies seen in the X-ray 
data reflect the local gravitational acceleration (assuming linear theory). 

We again note that the CG effect produces a dipole pattern on the sky of the form (see 
Eqn 6): 
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^ = (3 + q;)^cos(^) . (16) 
/ c 

Consequently the observed dipole will always be a coupling of the LSS and CG dipole 
anisotropies. 

The well known direct linear-theory relationship between the peculiar velocity and the 
absolute flux dipole is therefore (from Eqns 13, 14 & 15): 

(c.f. Boldt 1987, Lynden-Bell et al 1989, Miyaji 1994). To express the hnear velocity in 
terms of the LSS dipole anisotropy A we recall the Olbers integral for / (Treyer et al 1998, 
Eqn 11): 

Px{Q)rHQ{q,Zmax) ^^g^ 

47r 

where the Hubble radius rn = c/Hq and Q{q, Zmax) = ((1 + ^maxY '^'^ ~ ~ 2.5). Since 
the monopole Mj. — Anl and D^/M^ = A/3, for the fiducial values of q and z^ax used here 
we arrive at: 

Q0.6 

Ivl = 2.2 X lO^A-4- kms-^ . (19) 



7. Comparison with observed bulk motions 

As discussed above, we estimate that the bright sources we remove from the HEAO 
data are distributed to a distance ~ 60h^^Mpc, we can therefore compare our XRB dipole 
measurement with direct measurements of bulk flows over this scale. 

Several studies provide generally consistent estimates of the bulk flow of a ~ 50h~^Mpc 
radius sphere; 305 - 370(±110) kms^^ (Mill catalogue, POTENT, Dekel et al 1999), 
~ 300 kms"^ (SFI data. Dale et al 1999, GiovaneUi et al 1997), ~ 250kms-i (SNla data, 
Riess et al 1995). The directions of these flows are summarised in Figure 13. Using a 
crude mean of these numbers we estimate vqq ~ 300 ± 100 km s~^ The range of dipoles 
measured in the Total band (after removal of the dominant X-ray sources from within 
~ 60 h~^ Mpc) is 0.0023 < A < 0.0095 (depending on the measurement method used and 
corrections for the Galaxy and GG effect). This anisotropy is at most 2 times larger than 
the expected X-ray CG dipole. Applying Equation 19 the dipole measurements imply that 
1/7.1 < fio^/6a;(0) < 1/1.7. The favoured Method 2 dipole amplitude given in the last row 
of Table 2 is A = 0.0065 which yields Og V6^(0) = 1/4.8. 
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The quantity fl^^/bxiO) has also been estimated from studies of the X-ray selected 
AGN dipole under certain assumptions about local dynamics. Generally f2Q^/6^(0) ranges 
from 1/3.5 if all the local gravitational acceleration is assumed to arise from the volume 
with R ^ 45 h~^Mpc, to 1/7 if only half of the acceleration arises from within this volume 
( IMiyaji 199| ). Using the new IRAS-PSCz survey, Schmoldt et al (1999) predict that some 



65% of the LG acceleration is generated within 40h~^Mpc and that convergence is not 
reached until ~ 140h~^Mpc. Therefore Q^'^/bx{0) is almost certainly larger than 1/3.5 
using this method. These results are in good agreement with our above constraints from 
bulk flows and the HEAO dipole. The observed HEAO dipole therefore appears to be quite 
compatible with current measurements of the bulk flow of the local ~ 60h~^Mpc volume. 
The dominant population of X-ray emitters (AGN) in the 2-10 keV band is then highly 
biased with respect to other tracers, e.g. optical or IRAS galaxies. 

Over larger scales (~ 100 — 150h~^ Mpc) there is less agreement on the reality of bulk 
flow measurements. For example, the work of Lauer & Postman (1994) has suggested, with 
much controversy, that the Local Group has a motion relative to the z < 0.05 Abell cluster 
frame of 561±284 kms~^ in a direction I = 220°, b = — 28°(±27°). Assuming a dynamical 
origin of the observed CMB dipole this implies that the Abell cluster frame (to z = 0.05) 
is itse// moving in bulk with respect to the CMB frame with velociy of 689±178 kms^^ 
towards / = 343°, b = 52°(±23°). If correct this could imply that ~ 50% of the Local Group 
motion is due to matter on scales > lOOh^^Mpc. This specific result has been strongly 
refuted by several other works (e.g. Riess et al 1995, Giovanelli et al 1997, Hudson et al 
1999, Muller et al 1998, Dale et al 1999). More recently however, independent observational 
evidence for bulk motion over these scales has emerged, in the range of ~ 600 — 700 km s~^ 
(e.g. Hudson et al 1999, Willick 1998). All such studies however obtain directions for these 
motions greater than 60° away from the LP result, and are themselves highly fraught with 
potential systematic effects. 

In criticism of these results it can be noted that there is an inconsistency between (for 
example) the Lauer & Postman measurement and the results of gravity dipole estimation 
using galaxy catalogues. For example, the results of Strauss et al (1992) using the 1.2Jy 
IRAS redshift survey found an extraordinary convergence of the direction of the inferred 
gravity dipole out to ~ 20,000kms~^. This convergent dipole direction is only some 20° 
from the CMB dipole direction. (There are good arguments why the velocity dipole of the 
Local Group is not necessarily converged until 2 ~ 1 ( Pr'eacock 1992| ), but that does not 
preclude a genuine convergence in a smaller volume). Recently the IRAS-PSCz survey has 
largely confirmed these observations (Schmoldt et al 1999). 

Scaling our above estimates for the HEAO XRB dipole anisotropy we predict that, if all 
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X-ray sources within ~ 100 — 150h~^Mpc were removed then a 700 km s~^ bulk flow would 
correspond to 0.0054 < A < 0.0225 (assuming the measured range of allowed fl^'^ /bx:{0)). 
This would be some 1.5 — 5.0 times larger than the expected X-ray CG dipole amplitude. 

The directions of both the bulk flow estimated from other works, and our present dipole 
measurements are however scattered over a large area of the sky. Figure 13 summarizes 

most of these directions. As already mentioned, the LP flow direction is ^ 60° from all 
others, the Hudson et al (1999) result is also significantly further from the Solar CMB 
velocity direction. Interestingly the HEA01-A2 measurements appear to be somewhat 
intermediate to the LP result and the majority of the other, more local volume, estimates 
(SFl, Mill, PSCz). However, recalling that we estimate an XRB dipole direction error of at 
least ~ 30° for either method, then the HEAO dipole directions are not inconsistent with 
(for example) the SFI and Mill flows. 



8. Summary and Conclusions 

The HEA01-A2 X-ray data offers the best all-sky survey of baryonic matter to 
2; ~ 4 currently available. Although low level anisotropics in the X-ray sky background 
arise largely from Galactic contributions, relatively crude foreground removal clearly 
demonstrates the presence of extragalactic emission associated with well known large scale 
structure (Virgo, Coma, Centaurus/Great Attractor etc.) in the local universe. 

Quahtative comparison of the RASS 1.5 keV data with the 4 HEA01-A2 bands used 
here suggests that at least one third of the faint, unresolved HEA01-A2 structure may be 
Galactic in origin, and possibly associated with the Galactic Bulge. 

The local extragalactic hard X-ray emission is dominated by AGN and galaxy clusters. 
If we remove the flux associated with these sources to a flux hmit of 2 x 10~^^ erg s~^ cm~^ 
(2-10 keV) (which removes all sources more luminous than 5.2 x 10^^ erg s"^ (2-10 keV) 
out to 2; ~ 0.015) we measure a dipole anisotropy of A = 0.0023 to 0.0095 (depending 
on the method used and the details of the data processing). This range of anisotropy is 
consistent with our expectations (LPT 97) of comparable amplitude CG and LSS dipoles. 
It is significantly smaller than that measured in the RASS 1.5 keV all-sky data by Plionis 
& Georgantopoulos (1999). However, we have argued that the hard (> 2 keV) band XRB 
suffers less from Galactic contamination and have shown how removal of the foreground of 
bright sources reduces the dipole amplitude and shot noise (see also Treyer ct al 1998). 

We have derived the fully cosmological expressions for X-ray dipole anisotropy. Unlike 
the often used Euclidean case, the relationship of the local acceleration to the dipole 



-20- 



anisotropy is no longer straightforward. However, in the case of the current HEAO dataset 
we show that for reasonable choice of cosmology and matter density fluctuation power 
spectrum most of the dipole anisotropy arises from z ^ 0.1 and the low-redshift, linear 
theory approximation can be used. 

Using current estimates of the bulk flow of the local 60h~^ Mpc radius volume and our 
XRB dipole measurements we find that 1/7.1 < n° V^x(0) < 1/1.7. With our preferred 
dipole anisotropy measurement then QQ'^/bx{0) = 1/4.8. This implies that the population of 
X-ray sources is highly biased with respect to optical or infra-red selected objects. Studies 
of the dipole anisotropy of the local AGN distribution (Miyaji 1994), and other analyses of 
the HEAO data (Boughn et al 1998, assuming epoch-independent biasing) also yield high 
values. Interestingly our previous analysis of the angular power spectrum of the HEAO 
dataset (Treyer et al 1998) which included terms as high as Z = 20 yielded a present-epoch 
biasing factor 63.(0) ~ 1 — 2. The model fit to this data was however not particularly good, 
and the lower order harmonics (1 — 1 — 3) are better flt with a higher bx{0). If flo — 0.3 then 
the values of bx{0) estimated from the HEAO dipolc/bulk flows fall into this lower range, 
however our formalism is all based on an Q = 1, Einstcin-dc-Sitter cosmology. We also note 
that in all conventional models the bulk flow amplitude of a sphere with radius R drops 
with R (speciflcally, if P{k) oc /c" then Vbuik oc R~^^~^^^^^). Therefore, if we overestimate the 
volume within which we remove X-ray emission, but continue to apply the observed bulk 
flow amplitudes for a larger sphere, we will then underestimate flQ^/bx{0) from Equation 
15. 

If ~ 700 km s~^ bulk flows over 100 — 150 h^^ Mpc radius volumes did exist, as 
suggested by some studies (e.g. Lauer & Postman 1994 , Willick 1998, Hudson et al 1999), 
then we predict that an XRB dipole anisotropy of 0.0054 < A < 0.0225 would be seen after 
removing source emission within this volume. 

It is worth noting that we should not discount further complications such as a spatially 
varying local X-ray emmissivity to mass biasing. Indeed, in their study of the X-ray 
properties of the Great Attractor region and the Shapley supercluster Raychaudhury et al 
(1991) found that for these similarly massive regions the number counts of X-ray luminous 
clusters is quite different (Shapley having the most). This is suggestive of a spatial variation 
in cluster formation and brings into doubt the naive linear biasing scheme. 

To fully exploit this, or future, hard X-ray all-sky data for cosmological or dipole 
studies a better knowledge of the foreground contamination is essential. In particular, a 

signiflcantly more detailed model of the Galactic (and local, e.g. LMC, SMC) emission is 
needed. Probably the best way to achieve this will be through the use of softer band data (to 
provide spatial parameters) combined with point-by-point spectroscopic information which 
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will allow extrapolation to the harder, less contaminated X-ray bands. A spectroscopy 
capable all-sky imaging survey, such as that discussed by Jahoda (1998) would be well 
suited to this (see also discussion in Treyer et al (1998)). 

This research has made use of data obtained through the High Energy Astrophysics 
Science Archive Research Center (HEASARC), provided by the NASA Goddard Space 
Flight Center. CAS acknowledges the support of a National Research Council fellowship, 
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work was undertaken. CAS and MT thank the Racah Institute of Physics at the Hebrew 
University, Jerusalem, for its hospitality during completion of this work. 
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Dataset HEAO band Amplitude (A) Direction 



Galactic Mask 


Soft 


0.0132 


(325°, 51°) 




Jiard 


0.0084 


(355 ,45 ) 




Total 


0.0114 


(327°, 51°) 




R15 


0.0111 


(324°, 53°) 


Complete Mask 


Soft 


0.0036 


(343°, 35°) 




Hard 


0.0060 


(29°, 2°) 




Total 


0.0036 


(345°, 43°) 




R15 


0.0030 


(346°, 48°) 


Complete Mask 


Soft 


0.0033 


(336°, 40°) 


3% Iwan model 


Hard 


0.0032 


(32°, 2°) 


removed 


Total 


0.0034 


(338°, 47°) 




R15 


0.0029 


(337°, 54°) 


Complete Mask 


Soft 


0.0030 


(3°, -3°) 


CG effect 


Hard 


0.0048 


(41°, -27°) 


removed 


Total 


0.0027 


(8°, 5°) 




R15 


0.0022 


(15°, 0°) 


Complete Mask 


Soft 


0.0025 


(1°, -4°) 


CG effect & 


Hard 


0.0045 


(18°, -44°) 


3% Iwan 


Total 


0.0023 


(6°, 6°) 


removed 


R15 


0.0017 


(15°, 1°) 



Table 1: Spherical harmonic (Method 1) dipole measurements of HEA01-A2 data, directions 
in Galactic coordinates. 
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Dataset Amplitude (A) Direction 

Galactic Mask 0.0225 ± 0.0005 (335° ± 24, 46° ± 12) 

Complete Mask 0.0095 ± 0.0005 (344° ± 24, 35° ± 12) 

Complete Mask 0.0090 ± 0.0005 (327° ± 24, 25° ± 12) 

3% Iwan model removed 

Complete Mask 0.0085 ± 0.0005 (353° ± 24, 1° ± 12) 

CG effect removed 

Complete Mask 0.0065 ± 0.0005 (342° ± 24, 7° ± 12) 

CG effect & 

3% Iwan removed 

Table 2: Least squares fit (Method 2) dipole amplitudes of Total band HEAO data. Errors 
correspond to finite search grid scales. 
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Siimilatioii Dipole iiieasure Mean Amplitude (A) Mean sepii. from (2(35°. 18°) 

Complete Mask, Method 1 0.00135 ± 0.00031 

randomized fluxes Method 2 0.0022 ± 0.0014 



+ 4.2x10-3 CG effect Method 1 0.0024 ± 0.0007 31° ± 20° 

Method 2 0.0043 ± 0.0015 25° ± 15° 

Table 3: Flux dipole measurements from simulations. Amplitudes are taken as the mean 
over 10 realizations, errors correspond to 1 — cr standard deviations. 
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Fig. 1. — The full HEA01-A2 Total band data shown in Hammer- Aitoff Galactic coordinates 
projection. The Galactic Plane is clearly visible. Emission from the Large and Small 
Magellanic Clouds is visible towards the bottom of the plot, close to the southern ecliptic 
pole. 
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Fig. 2. — The mean counts/sec averaged over the time-wrapped Total band data (subject to 
Galactic masking) as a function of longitudinal pixel index. The least squares hnear model 
is shown as the solid line. 




Fig. 3. — The least squares linear trends to the time-wrapped HEAO data averaged at each 
longitudinal index using the complete source and Galaxy mask. Solid line - Soft band, dotted 
hne - Hard band, short dashed hne - Total band, long-dashed hne - R15 band. 
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Fig. 4. — Harmonic reconstruction to l^ax — 10 (angular resolution ~ 18°) of the raw total 
band HEA01-A2 data in Galactic coordinates. Note the complete domination by the Galaxy, 
in particular X-ray emission from the region of the Galactic centre. The first solid contour 
is at the mean intensity, subsequent soild contours are above the mean (dashed below) with 
separation of ~ 72% of the mean. The uppermost solid contour is at ~ 720% of the mean. 
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Fig. 5. — Harmonic reconstruction to Imax = 10 in Galactic coordinates (angular resolution 
~ 18°) of the 4 HEA01-A2 bands with a Galactic Plane cut of |6| < 22° and all resolved 
Galactic sources removed as well as the LMC nad SMC regions. Masked regions are filled 
at the mean intensity of the unmasked regions. The first solid contour is at the mean 
intensity, subsequent solid contours are above the mean intensity, dashed contours are 
below the mean. Contour separation is 0.04 counts/s in all bands. Associated structures 
can be identified: Virgo/Coma (at 6^80° and spanning the Northern Galactic cap), 
Centaurus/Great Attractor (at I ~ 315° and b ~ 30°) 
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Fig. 6. — Harmonic reconstruction to l^ax = 10 of HEAO data with the complete mask 
apphed. Contours arc as for Figure 5 but now separated by 0.013 counts/sec. Some of 
the features seen may be associated with known structures. For example, the feature at 
/ ~ 80°, b ~ —45° may be associated wth the Pegasus cluster/supercluster which extends to 
z ^ 0.04, the soft feature at I ~ 10°, b ~ —30° may be Galactic (see text) although the Pavo- 
Indus-Telescopium structure lies in the same direction. The feature at I ~ 210°, b ~ —30° 
lies in the direction of Abell 530 and 400, and the feature at I ~ 280°, b ~ 60° matches the 
direction of the Leo cluster /group. 
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Fig. 7. — In panels (a) and (b) the Imax = 10 harmonic reconstructions of example total 
band 'noise' skies are plotted for the two mask cases in Figures 5 & 6, using identical contour 
levels to these figures. The spherical harmonic coefficients are generated from Monte-Carlo 
realisations as described in Section 5 below, and reflect a random field drawn from the 
appropriate angular power spectrum in both cases (Section 4). Sky masking is also applied 
to these noise realisations. 




Fig. 8. — Harmonic reconstruction to = 10 of the RASS 1.5 keV data. The complete 
HEA01-A2 mask (Galactic and extragalactic sources removed) has been appUed to the data. 
Horizontal dotted lines at 6 = ±20° delimit the Galactic plane mask. Gontours are spaced 
at ~ 1.7% of the monopole, uppermost contour is at ~ 17% of the monopole 
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Fig. 9. — The normalized spherical harmonic power spectrum of the HEAOl XRB is plotted 
to / = 20 for Soft (dotted), Hard (short-dashed), Total (long-dashed) and R15 (dot-dashed) 
colours. The Galactic mask and extragalactic mask has been applied, corresponding to source 
removal to a flux limit of 2 x 10~^^ erg s~^ cm'^ (2-10 keV). The solid line shows a fiducial 
'best-fit' model for the Total band, as described in ( [Ireyer et al 19981) : low density CDM 
P(A;) with pure X-ray source density evolution and constant biasing (^^^(O) = 1), including 
the expected noise (see also §6). 
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Fig. 10. — The mean angular power spectrum over 10 realizations for each HEAO colour 
is plotted for 'noise' skys. Error bars correspond to the 1 — a scatter expected between 
individual realisations. 
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Fig. 11. — The Method 1 flux dipole amphtude in the Total band is plotted as a function 
of increasing Galactic correction (Iwan model, normalised to % of monopole). The upper 
curve (solid squares) corresponds to the dipole without removal of the expected CG effect, 
the lower curve (open squares) has had the expected CG effect removed following the Iwan 
model Galactic correction. 
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Fig. 12. — The growth of the rms flux dipole is plotted as a function of outer radius cut- 
off Zmax for a low density CDM P{k) with pure density evolution and constant biasing 
as described in ( [Treyer et al 1998]) . The solid curve assumes Zmin = and a flux cutoff 
S. 



cut 



2 X 10 erg s cm , while the thin lines simply assume a radial cutoff of z„ 



(upper dotted line) and 0.015 (lower dashed line) respectively (for both cases Wcut{z] 



= 0.01 
= 1. 
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Fig. 13. — A summary of key directions and dipole results plotted in equal area Galactic 
coordinates. Alphanumeric keys correspond to: LP - direction of Lauer & Postman (1994) 
bulk motion in CMB frame, CMB - direction of motion of the Sun in CMB frame (direction 
of expected CG effect) ( [Lineweaver et al 1996|) , SFI - bulk flow direction of 50h~^Mpc volume 
by Dale et al (1999) using SFI data (Giovanelli et al 1997), PSCz - predicted bulk flow in 
same volume from IRAS PSCz data (Schmoldt et al 1999), Mill - bulk flow in same volume 
from Mark III, POTENT catalogue (Dekel et al 1999), SNIa - bulk flow in same volume from 
Riess et al (1995), LPlOk - bulk flow of 150 h'^Mpc volume Tully-Fisher dataset (Willick 
1998), SMAC - bulk flow of 120h~^Mpc volume (Hudson et al 1999). Arrows correspond 
to Method 1 HEA01-A2 dipoles with increasing Iwan Galactic corrections (1-13%), upper 
(higher latitude) and lower arrows correspond to upper and lower curves in Figure 10 (without 
and with removal of the expected CG effect, respectively). Upper and lower asterix's label 
the corresponding Method 2 dipoles with 3% Iwan model corrections (Table 2). 



